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Abstract. The first part of the present paper is devoted to a systematic construction 
of continuous-time finite-dimensional integrable systems arising from the rational su(2) 
Gaudin model through certain contraction procedures. In the second part, we derive an 
explicit integrable Poisson map discretizing a particular Hamiltonian flow of the rational 
su(2) Gaudin model. Then, the contraction procedures enable us to construct explicit 
integrable discretizations of the continuous systems derived in the first part of the paper. 

1. Introduction 

The models introduced in 1976 by M. Gaudin [14] and carrying nowadays his name 
attracted considerable interest among theoretical and mathematical physicists, playing a 
distinguished role in the realm of integrable systems. 

The Gaudin models describe completely integrable classical and quantum long-range 
interacting spin chains. Originally the Gaudin model was formulated [14] as a spin model 
related to the Lie algebra 51(2). Later it was realized [15,20] that one can associate such 
a model with any semi-simple complex Lie algebra g and a solution of the corresponding 
classical Yang-Baxter equation [5,37]. Depending on the anisotropy of interaction, one 
distinguishes between XXX, XXZ and XYZ models. Corresponding Lax matrices turn out 
to depend on the spectral parameter through rational, trigonometric and elliptic functions, 
respectively. Both the classical and the quantum Gaudin models can be formulated within 
the r-matrix approach [34]: they admit a linear r-matrix structure, and can be seen as 
limiting cases of the integrable Heisenberg magnets [39], which admit a quadratic r-matrix 
structure. 

In the 80-es, the quantum rational Gaudin model was studied by Sklyanin [38] and 
Jurco [20] from the point of view of the quantum inverse scattering method. Precisely, 
Sklyanin studied the su(2) rational Gaudin models, diagonalizing the commuting Hamil- 
tonians by means of separation of variables and underlining the connection between his 
procedure and the functional Bethe Ansatz. In [12] the separation of variables in the 
rational Gaudin model was interpreted as a geometric Langlands correspondence. On the 
other hand, the algebraic structure encoded in the linear r-matrix algebra allowed Jurco 
to use the algebraic Bethe Ansatz to simultaneously diagonalize the set of commuting 
Hamiltonians in all cases when g is a generic classical Lie algebra. We have here to men- 
tion also the the work of Reyman and Semenov-Tian-Shansky [34]. Classical Hamiltonian 
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systems associated with Lax matrices of the Gaudin-type were widely studied by them in 
the context of a general group-theoretic approach. 

Some others relevants paper on the separability property of Gaudin models are [1,9, 
10,17,21,39]. In particular, the results in [9,12] are based on the interpretation of elliptic 
Gaudin models as conformal field theoretical models (Wess-Zumino-Witten models). As 
a matter of fact, elliptic Gaudin models played an important role in establishing the 
integrability of the Seiberg-Witten theory [36] and in the study of isomonodromic problems 
and Knizhnik-Zamolodchikov systems [11,30,35]. Important recent work on (classical and 
quantum) Gaudin models includes: 

• In [10] the bi-Hamiltonian formulation of sl{n) rational Gaudin models has been 
discussed. A pencil of Poisson brackets has been obtained that recursively defines 
a complete set of integrals of motion, alternative to the one associated with the 
standard Lax representation. The constructed integrals coincide, in the si{2) case, 
with the Hamiltonians of the bending flows in the moduli space of polygons in the 
euchdean space introduced in [22]. 

• In [18] an integrable time-discretization of su(2) rational Gaudin models has 
been proposed, based on the approach to Backlund transformations for flnite- 
dimensional integrable systems developed by Sklyanin and Kuznetsov [24]. 

• Integrable g-deformations of Gaudin models have been considered in [4] within the 
framework of coalgebras. Also the superalgebra extensions of the Gaudin systems 
have been worked out, see for instance [7, 13,29]. 

• The quantum eigenvalue problem for the gl(n) rational Gaudin model has been 
studied and a construction for the higher Hamiltonians has been proposed in [41]. 

• Recently a certain interest in Gaudin models arose in the theory of condensed 
matter physics. In fact, it has been noticed [2,33] that the BCS model, describing 
the superconductivity in metals, and the sl{2) Gaudin models are closely related. 

Finally, we mention the so-called algebraic extensions of Gaudin models, which has 
been studied in [26,27,31] with the help of a general and systematic reduction procedure 
based on Inonii-Wigner contractions. These extensions constitute also the subject of the 
present paper, with a slightly different derivation. Suitable algebraic and pole coalescence 
procedures performed on the Gaudin Lax matrices with N simple poles, provide various 
families of integrable models whose Lax matrices have higher order poles but share the 
linear r-matrix structure with the ancestor models. This technique can be applied for 
any simple Lie algebra Q and whatever the dependence (rational, trigonometric, elliptic) 
on the spectral parameter be. The models characterized by a single pole of increasing 
order and with q = su(2), will be called here the one-body su(2) tower. The base of 
the rational tower (corresponding to = 2) is nothing but the Lagrange top, a famous 
integrable system of classical mechanics. The many-body counterpart of the Lagrange 
top is called a Lagrange chain, it is a homogeneous integrable chain of Lagrange tops with 
a long-range interaction. On the other hand, the first element of the elliptic one-body 
su(2) tower is a particular case of the (three-dimensional) Clebsch system, describing the 
motion of a free rigid body in an ideal incompressible fluid, see [32]. 
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A systematic approach to algebraic extensions of Gaudin models appears independently 
in [8] and [26]. We remark that in [8] only sl{n) Gaudin models are considered and no 
r-matrix formulation is provided, as opposed to [26]. 

The present paper is devoted to the construction of an integrable time discretization 
of the rational su(2) Gaudin model and its one-body and many-body extensions. The 
theory of integrable maps got a boost when Veselov developed a theory of integrable 
Lagrangian correspondences [42], - symplectic multi- valued transformations possessing 
many independent integrals of motion in involution. Since then the theory of integrable 
discretizations has been substantially developed, a systematic presentation of the state 
of the art is given in [40]. Let us mention main common features of the discretizations 
found in the present paper: 

• They are genuine birational maps, not just correspondences. 

• They preserve an invariant Poisson structure but deform integrals, so that they are 
not Backlund transformations in the strict sense. However they can be interpreted 
as Backlund transformations for deformations of the original integrable systems. 

The paper is organized as follows. In Section [2] we recall the main features of the 
continuous-time rational su(2) Gaudin model in order to give a systematic construction 
of continuous-time one-body and many-body rational su(2) towers in Section El Section 
m is devoted to the explicit integrable time discretization of the rational su(2) Gaudin 
model. Then, in Section 0, suitable contraction procedures on the discrete Gaudin model 
allow us to provide integrable discrete-time versions of the whole one-body rational su(2) 
tower and of the Lagrange chain. In this context, the main goal is the derivation of 
continuous-time integrable systems and their discretizations: we say practically nothing 
about solving them. However, we always have in mind one of the motivations of integrable 
discretizations, namely the possibility of applying integrable Poisson maps for actual 
numerical computations. Finally, some concluding remarks are contained in Section [HI 

Let us present here our main results. Our departure point is the following Hamiltonian 
flow of the continuous-time rational su(2) Gaudin model: 

AiP + Ef=iy.,y.], l<^<N, (1) 

where G su(2), p G su(2) is a constant matrix, and pairwise distinct numbers Aj are 
parameters of the model. This flow admits N independent integrals in involution: 



N , > 

i^. = (p,y.) + ET^' i<^<^> 

— Afc — Ai 



(2) 



where (■, ■) denotes the scalar product in su(2) ~ M'^. 

An integrable explicit discretization of the flow ([T]) is given by 

y. = (i + eA.p)(i + eEf=iy.) y. (i + ^E7=iy.)"'(i + ^A.p)-\ (3) 

with 1 < z < iV. Here hat denotes the shift 1 1— > t + e in the discrete time eZ, where e is a 
(small) time step. The map ([3]) is Poisson w.r.t. the Lie-Poisson brackets on ©^su(2)* and 
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has independent and involutive integrals of motion assuring its complete integrability: 
H,{e) = ( p, y. ) + 5^ ^-^^ (^1 + ^A, A, ( p, p ) j - I 5^ ( p, [ y„ y,] ) . 

They are 0(e)-deformations of the original ones, given in Eq. ([2D. 

A contraction of simple poles to one pole of order A^ provides the integrable flow of 
the one-body rational su(2) tower, 

= [zo,Zi] + [p,Zi+i], 0<z<Ar-l, (4) 

with the convention zjv = 0. Its integrals of motion, 

^ k-l 

i/f ) = ( p, Zfc ) + - z„ Zfc_,_i ), < A; < AT - 1, (5) 

i=0 

are in involution w.r.t. the Lie-Poisson structure obtained through a (generalized) Inonii- 
Wigner contraction of ©^su(2)*, see eq. ( l20l) . An integrable discretization of the flow (jll) 
is given by the following map: 

N~i-1 

z, = (l + ezo)z,(l + £Zo)''+£[p,Zi+i]-2 ^ ad^p z,+„ {]<i<N-l.{Q) 

This map is explicit (one can compute Zj successively, from z = A^ — ltoi = 0), and 
Poisson w.r.t. the bracket (l20l) . preserving therefore the Casimir functions of this bracket. 
Additionally, it has A^ independent integrals of motion in involution, assuring its complete 
integrability: 

^ fc-l 2 

H'f\^) = ( P> Zfc ) + - ^( Zi, Zfe„i„i ) + - ( p, [zo, Zfc] ) + — ( p, p ) ^( Zi+i, Zfc_i ), 

1=0 i=0 

with < A; < A^ — 1 (these integrals are 0(e)-deformations of ([HI)). 

To stress the importance of the flow ([!]), we note that its simplest instance, correspond- 
ing to A^ = 2, describes the dynamics of the three-dimensional Lagrange top in the rest 
frame: 

Zo=[p,Zi], Zi = [zo,Zi], 

with zq G M? being the vector of kinetic momentum of the body, zi G being the 
vector pointing from the fixed point to the center of mass of the body, and p being the 
constant vector along the gravity field. The Lagrange top is a Hamiltonian system w.r.t. 
the Lie-Poisson bracket on e(3)*, with the Hamiltonian function 

f<y\ 1 

Its complete integrability is ensured by the second integral of motion h'^^ = ( p, zq ), and 
by the Casimir functions C^^ = (zo,zi) and C[^^ = |(zi,zi). The map 1^ for N = 2 
coincides with the integrable discretization of the Lagrange top found in [6]: 

zo = zo + £:[p,zi] , zi = (1 + £:zo)zi {l + ezo)'^, 
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with the deformed Hamiltonian function 

H?\£) = (p,zi) + ^(zo,zo) + |(p, [zo,zi]), 

(all other integrals remain non-deformed in this case). 

A contraction of = 2M simple poles to M double poles provides the integrable flow 
of the Lagrange chain, 



riii = [p, ail + 



^M 



l<i< M. 



Here (mj, a^) G e(3)* and /ij's are free parameters of the model. (In particular, for 
M = 1 and /ii = 0, one recovers again the Lagrange top, upon the re-naming zq ^ mi 
and zi ai.) The Lagrange chain possesses 2M independent integrals of motion in 
involution, given in Eqs. ( l42|43l) . An explicit discretization is given by 

rhi = {1 + e fiip) (l + eJ2jLiT^j^ (l + ^E^iimj) (1 + e/ii p)"^ + e [p, ] , 

Sii = {1 + e fiip) (l + eY^f^^mj^ ai (l + e Y,f=iT^j^ (l + £^yUiP)~\ 

with 1 < i < M. Expressions for the integrals of motion of this Poisson map are given in 
Eqs. ( 157|58|) . they are 0(e)-deformations of the integrals of the continuous system. 

2. The continuous-time rational su(2) Gaudin model 

The aim of this Section is to give a terse survey of the main features of the continuous- 
time rational su(2) Gaudin model. In particular, we give its Lax representation along with 
the interpretation of the latter in terms of the (linear) r-matrix structure. For further 
details we refer to [14,15,20,34]. 

Let us choose the following basis of the linear space 5u(2): 

-i \ 1 / -1 \ 1 / -i 

-i y) ' 

We recall that the correspondence 



= 2 I -i j ' = 2 U j ' = 2 V i 



1 



-1 — i — a? 



R-'^ 9 a = (a\ a^, a^) < — ^a=- .1,2 • s = a^da G su(2), 

^ ' ' ^ 2 \ — la^ + a^ la-^ J ^ " 

is an isomorphism between (su(2), [■,■]) and the Lie algebra {E?, x), where x stands for 
the vector product. (Here and below we assume the summation over the repeated Greek 
indices.) This allows us to identify vectors from with matrices from su(2). We supply 
su(2) with the scalar product ( ■, ■ ) induced from M'^, namely (a, b) = — 2tr(ab) = 
2tr(ba''^), Va, b G su(2). The matrix multiplication and the commutator in su(2) are 
related by the following formula: 

ab = -^(a,b)l + ^[a,b], Va,bGsu(2). (7) 

In particular, if ( a, b ) = 0, then ab + ba = 0. 

The above scalar product allows us to identify the dual space su(2)* with su(2), so 
that the coadjoint action of the algebra becomes the usual Lie bracket with minus, i.e. 
ad^ a = [ a, b ] = — adb a, with a, b G su(2) . 
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We will denote by {yf}^^^, I < i < N , the coordinate functions (in the basis cTq.) on 
the i-th copy of su(2)* in ©^su(2)*. So, = yfcTa ■ In these coordinates, the Lie-Poisson 
bracket on ©^su(2)* reads 

{y?,yj} = -kj^a(3^y7^ (8) 

with 1 < i,j < N. Here Sij is the standard Kronecker symbol and ea/^^y is the skew- 
symmetric tensor with €123 = 1- The bracket ([H]) possesses Casimir functions 

C^ = l{y^,y^), l<t<N. (9) 

Fixing their values, we get a symplectic leaf where the Lie-Poisson bracket is non-degenerate. 
It is a union of N two-dimensional spheres. 

The continuous-time rational su(2) Gaudin model is governed by the following rational 
Lax matrix from the loop algebra su(2)[ A, A^^]: 



N 



. ^ A - A, 

1=1 



where the Aj's, with Aj 7^ A^, 1 < i,k < N, are complex parameters of the model, and 
p G su(2) is a constant vector. This Lax matrix yields a completely integrable system on 
the Lie-Poisson manifold ©^su(2)*. In particular, its spectral invariants are in involution. 
This can be demonstrated with the help of a linear r-matrix formulation. We quote the 
following result [20]. 



Proposition 1. The Lax matrix [T^) satisfies the linear r-matrix relation 

{Cg{X)®l,l(^Cg{l2)}+[r{X- fi),Cg{X)(^l + l(^Cg{n)]=0, VA,/ieC, (11) 

with 

r(A) = -^(T« © (T„. (12) 
The r-matrix [W\) is equivalent to r(A) = —11/(2 A), where 11 is the permutation operator 

The spectral invariants of Cg{X) are the coefficients of its characteristic equation 
det(£g(A) — 1) = 0, which reads 



1 1 ^ 

y = 7(p.p) + 



4^^'^' 2 

i=l 



^ '' Hi ^ Ci 



[X-\ (A-A,)2J ■ 



Here Cj are the Casimir functions given in Eq. ([H]), whereas the functions 

N 

A,' — A,' 



N , > 



are the independent and involutive Hamiltonians of the rational su(2) Gaudin model. We 
shall focus our attention on Hamiltonians obtained as linear combinations of the integrals 
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Hi-. 



N ^ N N 



i=l 



2 ' ^ Xj — A, 



(14) 



i=l 



An important specialization of the Hamiltonian fll4p is obtained considering rii = A,, 
1 <i< N. It reads 

^ N N 

^G = ^$^(y.,y:,') + $^A,(p,y,). (15) 



i,j = l 
^+3 



From the physical point of view it describes an interaction of su(2) vectors (spins in 
the quantum case) with a homogeneous and constant external field p. One verifies by a 
direct computation that the Hamiltonian flow generated by the integral ([T^ is given by 



Eq. ( fT6i) admits the following Lax representation: 
with the matrix >Cg(A) given in Eq. (fTOj) and 



N 



H"^(A) = Ea^' M"^(A)=Ap + 5:y. 



1< i < A^. 



^e(A),-M[r^(A) 



(16) 



(17) 



3. Contractions of rational su(2) Gaudin models 

3.1. Contraction of the Lie-Poisson algebra ©^su(2)*. The following statement al- 
lows one to get the generalized Inonii-Wigner contraction of the direct sum of copies 
of su(2)* [19,26,43]. It shall enable us to construct the rational one-body su{2) tower in 
Subsection 13.31 See also [26, 27] for further details. 

Proposition 2. Consider the Lie-Poisson bracket ^ o/©^su(2)* ~ (M^)^ with coordi- 
nates {yj)jLi, and a linear map (R'^)^ (M^)^, (yj)jLi ^ {'^i)f=Q ! given by 

N 

z, = ^' J2^}yv 0<i<N-l, (19) 
i=i 

with pairwise distinct Vj G C and Q < d < 1 (contraction parameter). Then the bracket 
induced on (R^)^ with coordinates (zj)^g^ under the map 17^) is regular for •§ — >• 0, and 
tends in this limit to 

{z?,zf] = f*' ]X]t^^ (20) 

with < i, j < N — 1. We shall denote the Lie-Poisson algebra (d^) by Cn{su{2)*) . 
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Proof: Using Eqs. and (fT9i) we get: 

N 



a /9 



ri,m=l 



TV r ^7 

n=l 



The limit ?9 — leads to (|20|) . It is easy to check that the antisymmetric bracket fl20l) 
satisfies the Jacobi identity. 

□ 

The following N functions are Casimirs for the Lie-Poisson bracket ( l20l) : 

^ N-l 

C^f^ = -5^(zi,z^+fe_,_i), 0<fc<iV-l. (21) 



i=k 



We illustrate this construction by the cases of small N . For N = 2 the contracted 
bracket C2(su(2)*) reads 

' ^0 I ~ ~^a(3-f Zq, , I = —Sa/B-y , ^^Z^ , | =0. (22) 

This is the Lie-Poisson bracket of e(3)* = su(2)* (Bs Its Casimir functions are 

Cf = (zo,zi), Cf) = i(zi,zi). (23) 
For iV = 3 we get the contracted Lie-Poisson bracket C3(su(2)*): 

l^^O'^OJ' ~ I^^O'^f j" ~ ~^a/37 ^1 ) |^^0 ' ^2 ~ ~^a/37 ^2' (24a) 

'^r,^f} = -£«^^zj, |z-,4}=0, |4,;22^}=0. (24b) 

Its Casimir functions are 

= (Zo,Z2) + ^(Zi,Zi), = (Zi,Z2), = ^(Z2,Z2). 

The following result will be useful in the next Sections. 

Proposition 3. Let H,G be two involutive functions w.r.t. the Lie-Poisson brackets ^) 
on ©^su(2)*. IfH,G are the corresponding functions on Ci\f{su{2)*)) obtained from H,G 
by applying the map / fTPj) in the contraction limit d — > 0, then they are in involution w.r.t. 
the Lie-Poisson brackets (d^). 
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Proof: In the local coordinates {yf }a=i, 1 < « < A^, we have: 

. , ^ dH dG ( ^ ^\ ^dH dG . 



^ ^-1 fiH Fir 

1=1 n,m=0 " 

EdH dG ^ ^/ n\ 

n,m=0 " 
n+m<Af 

where the first term does not depend explicitly on the contraction parameter i). Perform- 
ing the limit — > we get {H, G} = 0. 

□ 

3.2. Contraction of the Lie-Poisson algebra ©^*^su(2)*. The following Proposition 
enables one to get a Lie-Poisson algebra given by the direct sum of M copies of Civ(su(2)*) 
directly from the Lie-Poisson algebra ©^*^su(2)* associated with a iVM-body Gaudin 
model. Its specialization to M = 1 is equivalent to Proposition [2l 

Proposition 4. Consider the Lie-Poisson brackets of ©^^^5u(2)* ^ (M^)^^^ with the 
coordinates {yj)fjf , and a linear map (R^)^*^ — (M^)^*^, (y^) ^-^ (zi^n), given by 

N 

Zi,n =^'Y. ^kn-l)+J yN(n-l)+j, 1 < U < M, < t < N - 1, (25) 



with pairwise distinct Uj G C and < 'd < 1. Then the bracket induced on (R ) with 
coordinates (zi^n) under the map l[25\) is regular for t? ^ 0, and tends in this limit t? — > 
to 

^'"'^■''"J \ i + 3>N, ^ ' 

with 0<i,j<N— 1 and 1 < n,m < M. We shall denote the Lie-Poisson algebra ^2^ 
by ©^C^(su(2)*). 

Proof: Using Eqs. ([H]) and (1231) we get: 

N 

|^i>5^fm|^ = ^N(n-l)+l^N(m-l)+k {^^(n-l)+^ ^Mm-l)+fc } ^ 

l,k=l 

N 
l,k=l 

N 

— Wn,m ta/37 ^ '^N{n-l)+l yN{n~l)+l ~ 

1=1 

1 + ] < N, 
0{^) i+3>N. 
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The limit ?9 ^ leads to pGl). 

□ 

The Lie-Poisson brackets ( l26l) have NM Casimir functions of the form ( 12T|) . 
A computation similar to the one in the proof of Proposition [3] leads to the following 
statement. 

Proposition 5. Let H,G be two involutive functions w.r.t. the Lie-Poisson brackets ^ 
on ©^*^5u(2)*. If H,G are the corresponding functions on ©*^CAr(su(2)*) obtained from 
H, G by applying the map in the contraction limit "i? — > 0, then they are in involution 
w.r.t. the Lie-Poisson bracket [2t 



3.3. The rational one-body su(2) tower. Our aim is now to apply the map ( fT9|) . in the 
contraction limit "i? — 0, to the Lax matrix ( ITOi) . in order to get a new rational Lax matrix 
governing the rational one-body 5u(2) tower. To do so a second ingredient is needed: as 
shown in [23,26] we have to consider the pole coalescence \i = -d vi, 1 < i < N . This pole 
fusion can be considered as the analytical counterpart of the algebraic one given by the 
map CT. 



Proposition 6. Consider the Lax matrix [W\) with \i = dvi, 1 < i < N . Under the map 
l[T^) and upon the limit — >• the Lax matrix [T^) tends to 



N-l 



MA) = P+ETiT' (27) 
while the Lax equation [T7\) turns into 



Cn{X)=[Cn{X),M^^\X) =- Cn{X),M'j;>{X) , (28) 

with 



N-l 



^n\X) = J2^^ M^j^\X) = Xp + zo. 
1=1 

The Lax matrix [2l\ ) satisfies the linear r-matrix relation fll\} with the same r-matrix 

m 



Proof: The first part of Proposition [6] can be proved by applying the map ( |T9l) and 
the pole coalescence Xi = 'dvi, 1 < i < A^, on Eqs. (fTOl) and (|T8l) . We get 

N iv iv_i j 

j=l J j=l i=0 ^ ^ 

and 

N „ N N-2 . „ X i+l Af-2 

- E - E E Ox «(^) — E ^ - x^-'(A). 



= 1 



j=l ■> j=l i=0 ^ ^ j=0 

N 



i=l 
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The fact that the Lax matrix (1271) satisfies the hnear r-matrix relation (fTTl) with the 
same r-matrix (fT2l) requires a longer but straightforward computation. We refer to [26, 
27, 31] for a detailed proof. 

□ 

The Hamiltonian flow described by the Lax equation fl2H|) is given by 



Zi = [zo,z,] + [p,Zi+i], 0<i<N~l, (29) 

with z^r = 0, while the characteristic equation of the Lax matrix det(£7v(A) — /i 1) = 
reads 

1 1 ^"^ H^^^ 1 ^"^ C^^^ 

~^ = 4 ( P' P ) + 2 X] 7^ 2 ^ ;\fc+Af+i ' 

A;=0 fc=0 

where the functions Cl^\ < k < N — 1, are the Casimir functions ( l2Ti) . while the 
functions 

^ fc-i 

i7f) = (p,z,) + -5^(z„z,_,_i), (30) 

are the independent involutive Hamiltonians of the rational one-body su(2) tower. 

Notice that it is possible to obtain the integrals ( l30l) using the map ( fT9l) . in the con- 
traction limit -(9 — 0, and the pole coalescence \i = dvi, 1 < i < N, from the integrals 
(|T3l). Let us flx i such that < i < - 1. We get 

N N -[ ^ i _ i 

k=l k=l j,k=l ^ 

N i-l N 



2 

A;=l m=Oj,fc=i 



(P>Zi) + ^E^ Z^, Zi_^_l ) = 



m=0 

In the above computation we have taken into account the polynomial identity 

i-l 

3 



m=0 

The contracted version of the Hamiltonian ( |T5l) is given by h[^\ namely the integral 
of motion generating the Hamiltonian flow given in Eq. (!29l) . while the contracted version 

of the linear integral Ylk=i = Ylk=ii P' y^ ) is given by H^^\ 

Let us remark that the involutivity of the spectral invariants of the Lax matrix Cn{X) is 
indeed ensured thanks to the r-matrix formulation ffTTj) . Their involutivity can be proved 
also without using the r-matrix approach, just by referring to Proposition [31 
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3.3.1. N = 2, the Lagrange top. Fixing = 2 in the formulae of the previous Subsection 
we recover the well-known dynamics of the three-dimensional Lagrange top described in 
the rest frame [3,6,16,23,34,40]. In other words the Lagrange top is the first element of 
the rational one-body su(2) tower. 

The Lagrange case of the rigid body motion around a fixed point in a homogeneous 
field is characterized by the following data: the inertia tensor is given by diag(l, 1, J3), 
I3 G M, which means that the body is rotationally symmetric w.r.t. the third coordinate 
axis, and the fixed point lies on the symmetry axis. 

The equations of motion (in the rest frame) are given by: 

Zo = [p,Zi], Zi = [zo,Zi], (31) 

where zq G is the vector of kinetic momentum of the body, zi G is the vector 
pointing from the fixed point to the center of mass of the body and p is the constant 
vector along the gravity field. An external observer is mainly interested in the motion of 
the symmetry axis of the top on the surface ( zi, zi )=constant. 

A remarkable feature of the equations of motion fl3T|) is that they do not depend 
explicitly on the anisotropy parameter 73 of the inertia tensor [6]. Moreover they are 
Hamiltonian equations w.r.t. the Lie-Poisson brackets on e(3)*, see Eq. f l22l) . 

The Hamiltonian function that generates the equations of motion flHTl) is given by 

i7f) = (p,zi) + ^(zo,zo), (32) 

and the complete integrability of the model is ensured by the second integral of motion 
Hq = (p, zo). These involutive Hamiltonians can be obtained using Eq. ( l30l) with 
N = 2, namely considering the spectral invariants of the Lax matrix £2 (A), see Eq. ( l27j) . 
The remaining two spectral invariants are given by the Casimir functions (l23ll . 

3.3.2. = 3, the first extension of the Lagrange top. Let us now consider the dynamical 
system governed by the Lax matrix flTTl) with = 3. The Lie-Poisson brackets are 
explicitly given in Eqs. f l24a|l24bl) . According to Eq. ( !30l) the involutive Hamiltonians 
are: 

Ho^ = {p,Zo), ijf^ = (p,zi ) + ^(zo,Zo), iff^ = (p,Z2) + (zo,Zi). 

Looking at the brackets fl24all24bl) and taking into account that zq and Z2 span re- 
spectively su(2)* and M^, we may interpret them as the total angular momentum of the 
system and the vector pointing from a fixed point (which we shall take as (0, 0, 0) G M^) 
to the centre of mass of a Lagrange top. Let us remark that zq does not coincide with 
the angular momentum of the top due to the presence of the vector zi. We think of zi, 
whose norm is not constant, as the position of the moving centre of mass of the system 

composed by the Lagrange top and a satellite, whose position is described by zi — Z2. 

(3) 

Here we are assuming that both bodies have unit masses. Notice that the integral HI 
formally coincides with the physical Hamiltonian of the Lagrange top flH^ where now the 
vector Zq is the angular momentum of system and the vector zi describes the motion of 
the total centre of mass. 
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According to Eq. fl^^ the Hamiltonian flow generated by the integral if} reads 

Zo=[p,Zi], Zi = [zo,Zi] + [p,Z2], Z2=[Zo,Z2]. 

We see that the vector zi does not rotate rigidly, though Z2 does. 

3.4. The rational many-body su(2) tower. The rational many-body su{2) tower may 
be constructed simply regarding the Lax matrix ( |271) as the local matrix of a chain of many, 
say M, copies of the Lie-Poisson structure Cn{su{2)*). Indeed the r-matrix formulation 
(JTTl) ensures that the Lax matrix 

with pairwise distinct poles fik of order describes an integrable system defined on 
(B^C]\i{su{2)*) with the same r-matrix formulation ffTTl) . See [26,31] for further details. 

Let us consider the special case N = 2, namely the Lie-Poisson algebra given by 
©*^c(3)*. The resulting integrable system has been called Lagrange chain in [26,28]. We 
now present a new derivation of such a system without using the r-matrix approach, but 
just considering the contraction procedure of a rational su(2) Gaudin model defined on 
©^*^su(2)*. According to Proposition H] the contraction of the direct sum of 2M copies of 
su(2)* (i.e. A^ = 2) leads to the Lie-Poisson brackets on ©^e(3)*. 

It is convenient to simplify the notation: 

zo,fc = nifc, zi,fc = afc, 1 < k < M. (33) 

We interpret = {ml, ml, ml) G and a^ = (a^,a^,a|) G as, respectively, the 
angular momentum and the vector pointing from the fixed point to the center of mass of 
the k-th top. The Lie-Poisson bracket on ©*^e(3)* is: 

|m^, m^l = -6k J e^p^ ml, |m^, a^j = -^j e^p^ al, |a^, a^j = 0, (34) 

with I < k,j < M. This bracket possesses 2M Casimir functions: 

gl') = (mfc,afc), Q(') = ^(afc,afc), 1 < k < M. (35) 

Using the notation introduced in Eq. (!33l) . the Lax matrix of the Lagrange chain reads 

M 

vCm,2(a) = p + 

i=l 

Let us now consider a rational su(2) Gaudin model with 2M poles. We have to apply 
the map defined in Eq. ( !25l) to the set of vectors {yi}f=i- 

(zo)i = mj = y2i + y2i_i, (zi)^ = a^ = {i^2i y2i + J^2i-i y2i-i), (37) 

with 1 < i < M. Moreover we define the following pole coalescence: 

>^2i = f^i + ^J^2i, \2i-1 = lii + i'2i-i, ^<i<M, (38) 

where the Aj are the 2M parameters of the rational su(2) Gaudin model. 



m, 



+ 



A - /ii (A - ^if 



(36) 
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Proposition 7. Consider the Lax equation [FT^ with the pole coalescence ^EEl- Under 
the map ( fTPI) and upon the limit i} it tends to 



with the matrix Cm,2{X) given by Eq. and 

AT 

A ci, 



M^m]2W = Y.-^ 

i=i 



/ij rrij 



Proof: We have: 



M 



Cg{\) 



+ 



M 



(39) 



(40) 



i=l 



y2t 



\X - fli - V2i~l X- fIi-^U2 



M 



p+E 



y2t-i + y2i 
A - /ii 



M 



E'&{l^2iy2i + 1^2i-iy2i-l) , ^/nN i?^0 ^ ... 
^^-^^ + o(^) — . /:m,2(a). 



A similar computation leads to the auxiliary matrices A^m2 ■'^l- (HO]) starting from 
the ones in Eq. ffTSj) . 

□ 



The Hamiltonian flow described by the Lax equation (139|) is given by 



m, 



_ P, ai ] + 



(41) 



with 1 < z < M, while the characteristic equation det(£A/^2(A) — /i 1) = reads 



/i 



1/ V 1^ 

t(p,p) + ^2^ 



k=l 



Rh 



+ 



+ 



(1) 



+ 



(2) 
k 



X-Hk (A-/ifc)2 (A-/ifc)3 (A-/ife) 



where the functions Q^,^"*, Qf^ are the Casimir functions (1351) . and the functions 



= ( p, mfc ) + 



1 



M 

(p,afe) + |(mfe,mfe) + ^ 



!afc,mj) (afc,aj) 



(42) 



(43) 



are the 2M independent and involutive Hamiltonians of the Lagrange chain. 

Notice that, as in the su(2) rational Gaudin model, there is a linear integral given 
by -Rfc = Sfcli(P;'^fc)- ^ possible choice for a physical Hamiltonian describing 

the dynamics of the model can be constructed considering a linear combination of the 
Hamiltonians Rk and Sk similar to the one considered in Eq. We have: 



(44) 



MM ^ M 

'Hm,2 = "^{fJ.kRk + Sk) = ^(p,yUfcmfc + afc ) + - y^(m^,mfc). 



k=l 



k=i 



i,k=l 
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It is easy to check that the integral (jSj) generates the Hamiltonian flow (|4T]) . If M = 1, 
the Hamiltonian (jUj) gives the sum of the two integrals of motion of the Lagrange top. 

We can construct the integrals of motion of the Lagrange chain also by using the 
Lie-Poisson map (1251) with the pole coalescence (l38ll directly in the Hamiltonians (JT3ll . 
according to 

Ri = lim [ H2i + H2i-i ] , Si = \im[i^ {v2i H2^ + z/2*-i H2i-i) ] . (45) 



4. Discrete- TIME rational su(2) Gaudin models 

The main goal of this Section is the construction of an integrable Poisson map dis- 
cretizing the Hamiltonian flow ffTUl) . We shall provide an explicit map approximating, for 
a small discrete-time step e, the time e shift along the trajectories of the equations of 
motion (fT6|) generated by the Hamiltonian function ( fT5|) . We have to remark that no Lax 
representation (hence no r-matrix formulation) has been found for this map. Its Poisson 
property and integrability will be proved by direct inspection. 

Proposition 8. The map 

: y..^% = {l + e A. p) (l + 5 Ef=iyj) yi (l + ^ E7=iy.) ~\l + e a, p)-^ , (46) 

with 1 < i < N and e ^ M., is Poisson w.r.t. the brackets ^ on ©^su(2)* and has N 
independent and involutive integrals of motion assuring its complete integrability: 

Hu{e) = ( P, y. ) + E (l + jA. A, ( p, p ) J - I 5^ ( p, [ y., y,] ) , (47) 

j^k j^k 

With l<k<N. 

Proof: Let us first notice that the map (l46l) reproduces at order e the continuous- 
time Hamiltonian flow (|T6ll . The map (H6l) is the composition of two non-commuting 
conjugations: = {T>^)2 o {V^)i, where 

■■ y^^y*={} + e Ef=,y,) y^{l + e Y^Uy) ' ' (48) 

(Df )2 : y* ^% = {l + e K P) y* {1 + eK P)-' , (49) 

with 1 < i < A^. Notice that {Vf)i o {V^)2 ^ {Vf)2 o (Pf )i. 

The Poisson property of the map V'^ is a consequence of the Poisson property of the 
maps {V^)i and {V^)2- In fact {V^)i is a shift along a Hamiltonian flow on ©^su(2)* 
w.r.t. the Hamiltonian J2f^k=iiyj^yk)- the other hand {pf)2 is a shift along a 
Hamiltonian flow on ©^su(2)* w.r.t. the Hamiltonian Y^=i{'Pi^ky*k) ■ Therefore the 
composition {'D^)2 o {'D^)i is a Poisson map w.r.t. the bracket ([8]). 

Let us now prove the complete integrability of the map fHUj) . We show that the func- 
tions fH7|) are indeed integrals of the map fHUl) . Their independence is clear, while their 
involution w.r.t. the brackets ([H|) is proved in Appendix 2. 
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Notice that the maps (HHl) . (I49l) imply, respectively, the following relations: 

N N 

{yhyj) = {yi,yj), yi + ^J^^yi^y^^^y^^ 2^^'^^''^''^' 

{P,yj) = {Pi,y*j), + [y„p] = y* + I Ai [p,y*] , (51) 

with 1 < i, j < N . The preservation of the functions ( H71) is demonstrated by the following 
computation: 

= (P'yn + E^^^(i + T^^^^-<P'P)j + |E<P'[y^'yl]) = 

N , \ / 2 \ ^ 

= (p,yO + E^)^^(i + T^^^^-(P'P)J-|E(P'[y^'y^-]) = ^^(^)' 

with 1 < k < N. Here we have used Eq. (|5Ti) in the first step and Eq. (|50i) in the second 
one. 

□ 

Using the discrete Hamiltonians (HTll we can compute the discrete-time version of the 
Hamiltonian ffTSl). It reads: 



1 ( \ 

Hgie) = ^ Xk Hk{e) = ^( p, y^ ) + - ^ ( yfc, y^ ) ( 1 + j (p, p ) j - 

AT 

~4 ^^^k-Xj) (p, [yfe,yi])- 

j,k=l 

Moreover we still have a linear integral given by Ylk=i^k{^) ~ X]feLi( P; yfe ); the 
continuous-time case. 



5. Contractions of discrete-time rational su(2) Gaudin models 

Performing the contraction procedures presented in Subsections 13. II and 13. 21 we can now 
construct the integrable discrete-time versions of the Hamiltonian flows fpUl) and fBT]) of 
the whole rational one-body su(2) tower and of the Lagrange chain. 
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5.1. The discrete-time one-body su(2) tower. The integrable Poisson map discretiz- 
ing the flow ( |29l) of the rational one-body su(2) tower is given in the following Proposition. 



Proposition 9. The map 

: z,- h-i^ z,- 

i=i 

with < i < N — 1 and e G M, zs Poisson w.r.t. the brackets IfW^) on Cn{su{2)*) and has 
N independent and involutive integrals of motion assuring its complete integrability: 

hI^\£) = {P,2k) + - ^{zi,Zk-i-i) + -(p, [zo,Zfc] ) + — (p,p) ^(zi+i,Zfc_i), (53) 

i=0 i=0 

with <k < N -1. 

Proof: Let us construct the map fl3^ through the usual contraction procedure and 
the pole coalescence Xi = ^ Ui, 1 < i < N, performed on the map ( H6l) . 

Consider the map (T>^)i in Eq. ( l48i) . Using the map ( fT9l) and assuming Aj = {^Ui, 
1 < i < N, we get: 

TV AT _^ 

k=l k=l 

= (1 + £:zo)zi (1 + £:zo)"\ 

with < 2 < — 1. Hence the contracted version of {'D^)i is given by 

Zi ^z* = (l + £Zo)z,(l + £Zo)-\ 0<i<N-l. 

On the other hand, a direct computation, with the help of Eq. (jTj), yields the contracted 
version of the map {'D^)2 in Eq. (H9l) : 

TV N 

Y^^'^fk = (l + £^z/fcp)y*(l + £^i/,p)-^ = 

k=l k=l 

N 



E E (-^y (1 + ^ ^ p) p' 



k=l j>0 

N-i-1 



< + 2 E (|)'ad^,z*^^. + 0(^). 



with < i < — 1. Performing the limit "i? we have: 

iV-j-l 

(Pf )2 : z* ^ z, = z* + 2 E (f )' ad^P 0<z<N-l. 

Now the composition {T>^)2 o {T^f)i is easily verifled to result in the map T)^ given in 
Eq. (1321) . The Poisson property of the map is a consequence of the one of the map 
in Eq. 



18 



MATTEO PETRERA AND YURI B. SURIS 



Next, we construct, by contraction of the functions (H71) . the integrals of the Poisson 
map ( |52l) . We know that fixing e = in Eq. (H7I) we recover the Hamiltonians (fT3l) of the 
continuous-time su(2) rational Gaudin model. Their contraction gives the Hamiltonians 
(1301) of the continuous-time rational one-body su(2) tower. Therefore it is enough to 
perform the contraction procedure just on the two e-dependent terms of the integrals 
(HZD. We have: 

N ^ i-1 

^^'^Hkie) = (p,Zi) + - ^(z„,z,„„_i) - 

k=l m=0 

N 

j,k=i 



1 e 

(P,Zi) + i:^{'2'm,'2'i-m-l) + - ( p, [ Zq, ] ) + 



m=0 
2 i-1 



m=0 j,fc=l 

with < i < - 1. 

The involutivity of the integrals {H'^^\e'))^~Q is ensured thanks to Proposition O 

□ 

Let us remark that the specialization to = 2 of the map fl52l) gives the integrable 
time-discretization of the Lagrange top found by A.I. Bobenko and Yu.B. Suris in [6]. 
According to Eq. ( |52l) it reads: 



zo = zo + £ [p,zi] , zi = (1 + £:zo)zi (1 -hezo) ^ (54) 

The above explicit map approximates, for small £, the time e shift along the trajectories 
of the Hamiltonian flow flHTj) . This distinguish the situation from the map in [25], where 
Lagrangian equations led to correspondences rather than to maps. 

The map ( 15^ is Poisson w.r.t. the bracket ( l22l) on e(3)* and its complete integrability 
is ensured by the integrals of motion 

iff = ( p, zo ) , iif ) (e) = i ( zo, zo ) + ( p, zi ) + I ( P, [ zo, zi ] ) . (55) 



A remarkable feature of the map flMj) is that it admits a Lax representation and the 
same linear r-matrix bracket ffTTj) as in the continuous case, see [6] for further details. 
The Lax matrix of the map is a deformation of the Lax matrix of the Lagrange top. 
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5.2. The discrete-time Lagrange chain. The integrable Poisson map discretizing the 
flow ( HTl) of the Lagrange chain is given in the following Proposition. 



Proposition 10. The map 

uii = (1 + e Hip) (^1 + 6 Y.jLi'^j 



m. 



(56a) 
(56b) 

with 1 < i < M and e eR, is Poisson w.r.t. the brackets [34^ on ©*^e(3)* and has 2M 
independent and involutive integrals of motion assuring its complete integrability: 



1 + e fiip) (l + eY^^i-^mj) a^ (l + e^f^^mj) (l + e/i^p) ^ 



Rk{e) 



(p,mfc)-^(p, mk,Y.jLi^j ) 



M r 



+ 



S,(e) 



(mfc,aj) 



1 + J /^fc /^j ( P, P ) 



1 + J /^fc ( P. P ) 



(mj,afc) 



1 + J ( P, P ) 



(57) 



(p,afc) + ^(mfc,mfc) ( l + ^/i^(p,p) ) -^(p, &k,Ylf=i™-j ) 



M 



/ifc /Xj ( p, p ) + 



/^fc ( P, P ) 



(5^ 



1< k < M. 



Proof: Using the map fl37p and the pole coalescence (1381) in the map 
we immediately obtain the contracted version of ("D^^^)!. It reads 



* = y*2i + y*2^-l = ( 1 + T.f=i^j ) mi ( 1 + ^ E 



a* = ^ {y2i y*2i + 1^21-1 y*2i^i) = [i + e Y.f=i™^j) ^ + ^ E^lii 

with \ <i < M . The same procedure leads to the contracted version of {T>1^^)2- It reads 

mi = y2i + y2i-i = (1 + e/iip) m* (1 + e/iip)""^ + 

+ e [p, (l + e/i,p) a* (l + £^,p)~^] +0(?9), (60a) 

Sii = ?9(z^2iy2i + «^2i-iy2i-i) = (1 + e/iip) a* {1 + e iupY^ + Oi^d). (60b) 

Performing the limit ?9 — * in Eqs. f l60a|60bl) and combining the resulting equations with 
the maps in Eqs. (I59all59bp we obtain the map (I56all56bl) . Its Poisson property is ensured 
thanks to the Poisson property of the map fHUl) . 

The construction of the discrete Hamiltonians (I57|58p is similar to the one done for the 
continuous-time Lagrange chain. They can be obtained through the following formulae 



, m 



with N = 2M 

(59a) 
(59b) 
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by a straightforward computation: 

Ri{e) = \im[H2i{e) + H2i-i{e)' 



Si{e) = hm ['&{u2i H2i{e) + V2i~i H2i-i{e)) ], 

being the Hamiltonians fl47D . 
Let us finally notice that the Hamiltonians fl57ll58p are in involution w.r.t. the brackets 
thanks to Proposition 

□ 

The discrete-time version of the Hamiltonian ( l44l) is given by 

M 



k=l 

M 



1 ( \ 

^( p, /ife rrifc + afe ) + - ^ ( mj, mfc ) M + — ^jTXfe ( p, p ) j - 

k=\ j,k=l ^ ^ 

M 

-J^(^^k-^ij) ( P, [ mfc, nij ] ) - - ( p, Ej'ilafc, Ej^imj ) + 



j,fc=i 



A/ 



A/ 



+ y ( p, p ) XI ^ + T ^ P' p ^ 5Z ^ ^• 



j,fc=l 



j,k=l 



Notice that we still have the linear integral Ylk=i ^ki^) = J2k=i{ P' ™fc )• 



6. Concluding remarks 

We presented a systematic construction of finite-dimensional integrable systems shar- 
ing the same linear r-matrix bracket with the rational su(2) Gaudin model. The resulting 
one-body and many-body integrable systems are obtained through suitable algebraic con- 
tractions of the Lie-Poisson structure of the ancestor model. We called these families of 
integrable systems su(2) towers. The three-dimensional Lagrange top is the first element 
of the rational one-body 5u(2) tower. The many-body counterpart of the Lagrange top, 
called Lagrange chain, is also presented and its Lax representation is given. 

In the second part of the paper we derived an explicit integrable Poisson map dis- 
cretizing a Hamiltonian fiow of the rational su(2) Gaudin model, thus providing a new 
integrable discretization of such a model. Then, the contraction procedures enable us to 
construct integrable discrete-time versions of the of the rational su(2) tower and of the 
Lagrange chain. 

The main open problem connected with this work is to find Lax representations (and 
then their r-matrix interpretation) for all the integrable Poisson maps introduced here 
(actually the only case for which the Lax representation is known is the discrete-time La- 
grange top considered in [6]). These structures will allow to avoid a brute force verification 
of the integrability, which we had to perform here. Of course, finding a Lax representation 
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for the discrete-time rational 5u(2) Gaudin model would yield the corresponding results 
for all the contracted systems. 

Also the following problem deserves further investigations. It is well-known that the 
continuous-time rational Gaudin models, as well as the one-body and many-body towers 
[31] described in the present work, admit a multi-Hamiltonian formulation [10]. Finding 
a multi-Hamiltonian formulation of our discrete-time maps is an open challenge. 
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Appendix 1: Visualization 

As shown in Proposition [9] the integrable Poisson map ( l52l) discretizing the rational 
one-body su(2) tower is well defined and can be easily iterated. We present here its 
visualization in the case N = 2 (i.e. the Lagrange top) and = 3 (i.e. the first extension 
of the Lagrange top). 

The input parameters are: the intensity of the external field, p; the discretization 
parameter, e; the number of iteration of the map, J\f; the initial values of the coordinate 
functions: (zo(0), Zi(0)) for A^ = 2 and (zo(0), Zi(0), Z2(0)) for A^ = 3. 




The first plots refer to the case N = 2. The output is a 3D plot of A/" consequent 
points [zl, zf, zf), describing the evolution of the axis of symmetry of the top on the 
surface (zi,zi )=constant. These plots show the typical (discrete-time) precession of the 
axis. 

The second ones refer to the case A^ = 3. The output is a 3D plot of A/" consequent 
points {z2, Z2, Z2), describing the evolution of the axis of symmetry of the top on the 
surface ( Z2, Z2 )=constant and Af consequent points {zl — z\,z1 — z^, z\ — ) describing 
the evolution of the satellite. 
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Let us also give a visualization, for M = 2, of the integrable discrete-time evolution of 




The input parameters are: the intensity of the external field, p; the values of the 
parameters fii and fi2', the discretization parameter, e; the number of iteration of the 
map, Af; the initial values of the coordinate functions, (mi(0), ai(0)) and (1112(0), a2(0)). 

The output is a 3D plot of Af + Af consequent points {a\,a1,al) and {al,a2,al) de- 
scribing the evolution of the axes of symmetry of the tops respectively on the surfaces 
( ai, ai )=constant and ( a2, a2 )=constant. 

Appendix 2: Proof of the involutivity of the functions {Hk{e)}^^^ 
Let us write the functions {Hk{e)}^^^ given in Eq. (H7j) in the following way: 

Hk{e) = K-'-hl + J {p,p) hi l<k<N, 
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where 



N 



hi 



hi 



{yk,yj) 



j=i ^'^ ~ 



N 



hi = X]^P'[y^''yj])' 



N 



(61a) 



(61b) 



(61c) 



In the following computations we shall use the Lie-Poisson brackets ([H]). We have: 

{H,{e),H,{e)] = [hlX] ' ^ {[hlM] + [hlX]) + 

+ j [(P,P) {[hlM] + [hlX]) + [hlM]] - 

- J {{hi h^} + {hi, h]}) + 1^ ( p, p h^} . (62) 

We already know that {h^.^h'^} = 0, 1 < k,i < N, since the integrals {h^}^^^ are the 
ones of the continuous-time su(2) rational Gaudin model. Let us compute the remaining 
brackets in Eq. ( |62l) using Eqs. fl61ay61by61cl ) and assuming k ^ i. Notice that in the 
brackets h}} + {hi, /i^} and {h^., hf} + {hi, /i^} we shall exphcitly write the order of 
IpI = ( P) P )^''^ = P appearing in the computation. 
At order e we have: 



[{hlhl} + {hl,h^}]^^^^^. 



N N 



Afc - Aj 



{ykVj.yfyi} 



A,; — Aj 



j=l 1=1 

j^k l^ti 

^ 1 

-P^ ei3pa eaa^ X - X- ^^^^ ^ " 



{y^yj^vry?} 



j=l " 

N 



N 



eppa e<xa7 X- - A- ^ ■ 



24 MATTEO PETRERA AND YURI B. SURIS 

The above expression vanishes if we swap the indices a and 7 in each second term in 
three brackets. Then we have: 

N 

1=1 

that vanishes due to the properties of the tensor e^/j-y. 
At order we get: 



[{ft2,A|} + {ALA°}]o,|p|,= 



1=1 j=i 



= -P ea/37 T ^ViVk + -^ ^ViVk 

that vanishes swapping the indices 7 and (3 in the second term. Moreover, 



N N \ \ I \ \ 

(A.-A,)(A,-A,)(A,-A,) 

TV 



On the other hand: 



iV TV AT 



{hl,hl} = ead^e ^py. Y^^ivlv]. Vi Vi] = P" P" ea/37 vl Vi, 



j=i 1=1 j=i 

j^k l^t j^k 



where we have used the properties of the tensor Safj^y- Hence we get: 



(p,p) m,h^}+{hih^})+{hihi}=o. 
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At order we have: 

{hlh^} + {hlhl}^ 



- -/e,,.EE {ytyf,yfyi} + {y'^y^yryr} 

j=i /=i L ^ I I 

^ \ \ 



eppa eaaj 2^ T ^ {Vk yf yt + Z/J ) " 

^ \ \ 

- e,,. 6,., ^ (y,^ + y7 yt) ■ 

- P' 6,,. 6<,„, J] {yi y] yt + y] y^ ^ -"^ 



3^k 

The above expession vanishes if we swap the indices a and 7 in each second term in the 
three brackets. 

Finally, at order e^, we get: 

^ ^ >^k^3^i ^ hXjXf XkX]Xi 



3¥^k 



{Xi — Xk){Xk — Xj) {Xi — Xj){Xk — Xi) {Xi — Xj){Xk — Xj) 



A direct computation shows that the expression in the square brackets vanishes. 
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